Critical behavior of the Ising model in annealed scale-free networks 
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We study the critical behavior of the Ising model in annealed scale-free (SF) networks of finite 
system size with forced upper cutoff in degree. By mapping the model onto the weighted fully 
connected Ising model, we derive analytic results for the finite-size scaling (FSS) near the phase 
transition, characterized by the cutoff-dependent two-parameter scaling with four distinct scaling 
regimes, in highly heterogeneous networks. These results are essentially the same as those found for 
the nonequilibrium contact process in annealed SF networks, except for an additional complication 
due to the trivial critical point shift in finite systems. The discrepancy of the FSS theories between 
annealed and quenched SF networks still remains in the equilibrium Ising model, like some other 
nonequilibrium models. All of our analytic results are confirmed reasonably well by numerical 
simulations. 
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I. INTRODUCTION 

Many aspects of our real world have been understood 
in the context of complex networks [l], Q and simple phys- 
ical models of critical phenomena on networks. Contrary 
to regular lattices in the Euclidean space, complex net- 
works are characterized by a highly heterogeneous struc- 
ture as manifested in broad degree distributions. Recent 
studies on equilibrium or nonequilibrium systems have 
revealed that the heterogeneity is one of essential ingre- 
dients determining the universal feature of phase transi- 
tions and critical phenomena [3j. 

The concept of the phase transition is well defined 
only in the thermodynamic limit where the system size 
is taken to infinity. So it is important to understand 
how finite-size effects come into play near the transition. 
Such a task for physical models on regular lattices has 
been successfully accomplished by the standard finite- 
size scaling (FSS) theory 0], based on the ansatz that 
a single characteristic length scale (correlation length) 
£ competes with the system's linear size L. Then, any 
physical observable depends only on a dimensionless vari- 
able £ = L/£ in the scaling limit. Near a second-order 
continuous transition, the correlation length diverges as 
£ ~ l e l" with the reduced coupling constant e and the 
finite-size effects become prominent. 

The FSS theory for complex networks can be formu- 
lated in a similar way: Since the Euclidean distance is 
undefined in complex networks, one may take the vol- 
ume scaling variable as l v — N/£ v with the system size 
N (the total number of nodes) and the correlated vol- 
ume £„. The correlated volume diverges £„ ~ near 
the transition (y = vd in d dimensional lattices). For 
example, the magnetization of the Ising model scales as 



m{e,N) = N-P^tpieN 1 ^), 



(1) 



where the scaling function ip(x) ~ 0(1) for small x and 



x@ for large x with the order parameter exponent /3. 

The FSS theory with a single characteristic size has 
been tested numerically in many systems (see Ref. Q and 
references therein). In particular, the exact values for the 
FSS exponent v are conjectured || by estimating the 
correlated volume (droplet) size for the nonequilibrium 
contact process (CP) and the equilibrium Ising model in 
random uncorrelated networks with static links, which 
are denoted as quenched networks. 

However, considering a highly heterogeneous scale- 
free (SF) network, one should take into account not only 
a broad degree distribution of P(k) ~ k~ x but also the 
upper cutoff k c in degree, which scales as k c ~ N 1 ^ . 
Without any constraint, k c is bounded naturally with 
w na t = A — 1. In general, one may impose a forced cut- 
off with uj > Wnat- In the thermodynamic limit, both TV 
and k c diverge simultaneously and u) sets a route to the 
limit. Therefore, one can suspect that the FSS theory 
may depend on the routes or equivalently on the value 
of to, especially for networks with a broader distribution 
for small A. 

For the quenched SF networks, it has been suggested 
that the FSS does not vary with to for a weak forced cut- 
off (to < A), which was confirmed numerically in various 
types of SF networks However, in the annealed net- 

works where links are not fixed but fluctuate randomly in 
time, it was rigorously shown that the CP model exhibits 
an anomalous FSS for any forced cutoff with 2 < A < 3 
where a hete roge neity (A) -dependent critical scaling ap- 
pears @, d, B EH • Moreover, the anomalous FSS is char- 
acterized by a cutoff(w)-dependent and two-parameter 
scaling with four distinct scaling regimes [Io| , in contrast 
to the cutoff-independent and single-parameter scaling 
with three scaling regimes in the standard FSS theory. 

The anomalous FSS of the CP in the annealed SF net- 
works gives rise to a natural question: What is the main 
ingredient causing the anomaly? Some possible guesses 
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may be a nonequilibrium feature of the CP, absorbing na- 
ture (vanishing activity) at criticality, or heterogeneity of 
networks [1, In this paper, we answer to this ques- 
tion by studying the Ising model, a prototype equilibrium 
phase transition model, in annealed SF networks. We 
find the same type of the anomalous FSS scaling (cutoff- 
dependent and two-parameter scaling with four distinct 
scaling regimes) for any forced cutoff with 3 < A < 5 
where the A-dependent critical scaling appears in the 
thermodynamic limit for the Ising version. In addition, 
the trivial shift of the critical point in finite systems adds 
one more complication on the critical FSS, though it does 
not cause any fundamental change. In summary, our re- 
sults may draw a general conclusion that the anomalous 
FSS scaling should appear in any critical system in the 
annealed SF networks for any forced cutoff (u> > w na t) 
with the degree exponent A such that a A-dependent new 
singularity arises in the physical quantities as N — ► oo. 

This paper is organized as follows. We define the Ising 
model on an annealed network in Sec. [IT] and show that 
it is equivalent to the Ising model on the weighted fully 
connected network. In Sec. IIII1 the FSS theory is devel- 
oped in various networks including SF networks, which 
is numerically tested in Sec. IIVI In Sec. [V] some effects 
of the sampling disorder are discussed. We conclude this 
paper with summary and discussion in Sec. IVII 



II. ISING MODEL ON ANNEALED NETWORKS 

An annealed network Gn is defined as an ensemble of 
all networks consisting of TV nodes which are assigned 
to a given degree sequence {k\, . . . , kjy}. An instance 
g € Gn is constructed by assigning hi stubs to each 
node i (1, . . . , N) and then completing edges by pairing 
the stubs randomly as in the uncorrelated configuration 
model HI 03. 

A network configuration g is conveniently represented 
by an adjacency matrix A(g) whose element Ay takes 
either 1 or if there is an edge between nodes i and j 
or not, respectively. In the ensemble Gn, the connecting 
probability Pijto find an edge between two nodes i and 
j is given by % Q]| 



Pij 



(2) 



with the mean degree z\ = ki/N . This expansion is 
valid when -C 1 for all i and j. 

The ferromagnetic Ising model on the annealed net- 
work Gn is defined by the Hamiltonian 



H[{s},g] = -jJ2 A ^(9)stSj - £/k 



(3) 



i<3 



where J > is a ferromagnetic coupling constant, Sj € 
{ — 1,1} is an Ising spin variable at node i, and hi is a 
local field at node i. In comparison to the model on 
a quenched network, a network configuration g is also 



fluctuating within Gn as well as the Ising spins. Ther- 
modynamic properties of the model is obtained from the 
partition function 



z = £ E ex p 

se6w {si} 



Ky^ j Aij(g)s i Sj + hjSt 



i<j 



> (4) 



where K = [3J and hi = (3hi with the inverse tempera- 
ture (3 = l/k B T. 

In terms of the connection probability in Eq. (J2), one 
can easily perform the partial summation over g to obtain 
that 

z = e n K 1 -pu) +pne Ksisi ] ip is< - ( 5 ) 

{si}i<j i 

Utilizing the identity e Ks = cosh A"+ssinh K for s = ±1, 
we find that 



Z = Z CX P ( E Qii s i s 3 + E ^ iSi 

{Si} 



(6) 



where Zq is an overall constant factor (not depending on 
{*}), 



and 



tanh Qi 



1 — + p^ cosh K 
coshQi ? - 



Pij 



sinh K 



1 - Pij + Pi 3 cosh K 



(7) 



As Qij is nonzero for any pair of the expression in 

Eq. © corresponds to the partition function of the Ising 
model on the fully connected network with the heteroge- 
neous coupling constants Qij. 

As p.^ — kikjj (Nzi) <C 1 for large N [l4j], one can ap- 
proximate Qij ~ Kkikj /(Nzi) with K = sinh A". Hence, 
in this paper, we focus on studying the Ising model on 
the fully connected network with the Hamiltonian Hf as 



i<3 



(8) 



This Hamiltonian was studied as a MF or annealed ap- 
proximation for the Ising model on quenched networks in 
the thermodynamic limit 0, ITU [T(| . 

For convenience, we rewrite Hf in a completed square 
form as 



pH f = - 



K 
2N Zl 



E ^ iS{ ) E 



J2hsi (9) 



and define the magnetic order parameter as 

M = J2hsi, (10) 



3 



with the order parameter density to = M/(Nzi), which 
is first suggested in [l7j and recently for both equilibrium 
and nonequilibrium models in [l8| . 

Now we derive the free energy as a function of M, 
which allows us to calculate thermodynamic properties 
not only in the thermodynamic limit by minimizing it 
with respect to M but also for finite size N, at least up to 
the leading order. After dropping the additive constant 
term in Eq. ([9]), the partition function, up to a constant, 
can be written as 



Z = 



E 



exp 



K 
2N Zl 



{Si} V 



- J " M 

where 
F(M,u) = - 



io ° du ~ ~ 

— exp[-F(M,u) 

-inn ^TTt 



K 



(11) 



2Nz x 



M 2 + uM - Mn[2 cosh(ufc; + hi)], 



where (), = -k YliQi denotes the average over nodes. In 
obtaining Eq. pip , we used the integral representation 
of the delta function 5{M) = f |e'" l] and the analytic 
continuation v = iu. 

The integration over u can be evaluated us- 
ing the steepest descent method, which yields that 
the free-energy function F(M) defined by Z = 
J dM exp(— F(M)) is given as 



F(M)~F(M,u ) + -In 



2n\F" 



(12) 



where F" — —Nkfsech 2 (uoh + hi) is the partial second 
derivative of F(M,u) with respect to u at Uq. The con- 
dition that the first derivative F'\ u — Uo = determines 
u = u a (M, {hi}) by 



M = Nki tanh(w fci 4 
or equivalently, uq — uo(rh, {hi}) by 



m = — ki tanh(wofci + hA 

Zl 



(13) 



(14) 



We remark that the second and high-order terms on the 
right hand side of Eq. (fT2")) can be neglected because they 
increase with system size N only logarithmically, in con- 
trast to the first bulk term. Moreover, the finite-size 
corrections near the transition are stronger than the con- 
tributions from these terms. 

It is convenient to use the free energy density function 
/(m) = F(M)/N which is 



f(m) 



ziK _ 



rh +z\Uofh — ln[2 cosh(itofcj + hi)]. (15) 



Then, the ensemble-averaged value (rh) can be calcu- 
lated, in the thermodynamic limit, as 



(rh) = — I dm me 
Z 



-Nf(rh) 



TOO, 



(16) 



where Z = Z/(Nz x ) = J dm e - Nf(rh) and m is the 
minimum point of /(to). Here, the higher-order finite- 
size corrections are again at most logarithmic. 

The spin magnetization rrii at node i can be obtained 
by differentiating the partition function in Eq. (1111) by 
the local field hi, which result in 



tanh(ufci + hi) 
tanh u ki + hi , 



(17) 



with u = u (m , {hi}). 



III. FSS THEORY IN ANNEALED NETWORKS 

We are now ready to investigate the bulk critical scal- 
ing and also the FSS of the Ising model on annealed net- 
works. First, we consider the simplest case of exponential 
degree distributions such as the Poisson distribution of 
the random network. Then, we proceed to discuss for the 
SF degree distributions with P(k) ~ fc~ A with an upper 
cutoff k r - N 1 ^. 



A. Exponential networks 

Consider exponentially bounded degree distributions 
such that the degree moments z n = kf are bounded for 
all n. The Poisson distribution for the random network 
and the Kronecker 5-function distribution [P(k) = 5^ 2 ] 
for the random z-regular network fall into this category. 

Taking the uniform magnetic field hi = h and expand- 
ing Eq. p4p for small uq and h, we get 

m = h + —u -^u 3 + o(ul,h 2 ,hul). (18) 
Then, the free energy density is given by 



a 2 b 4 



f (to) = — In 2 — ohm em" H m 

M 1 2 12 



where 



2 / i 4/4 

a = z 1 /Z2, o = z 1 Z4/z 2 



(19) 



(20) 



and the reduced inverse temperature e = [K — K c )/K c 
with the critical point 



K c = zi/z 2 . 



(21) 



Note that a = b = 1 and K c = l/z for the random 
z-regular networks. 
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At h — 0, the order parameter scales for e > as 



v/30/6 



(22) 



with the order parameter exponent (3 = 1/2. It is 
straightforward to derive the zero-field susceptibility x = 
d(m)/dh\~ h=Q ~ (2e)-Tfore > and x ^ (l-2/7r)(-e)-T 
for e < with 7 = 1. The average magnetization 
to = Wi is related to (to) through Eq. Q17p , which yields 
to ~ a(rh). 

With the free energy function given in Eq. (|19p. one 
can develop the FSS theory analytically. The full scaling 
functions for (to) and \ are derived in the Appendix. We 
only summarize the results below. The FSS form for the 
order parameter is given by 



(m(e,N)) = N'^^eN 1 ^;^^, 



(23) 



where = 1/2, the FSS exponent v = 2, and the scaling 
function ip is given by Eq. (|A3|) . The function arguments 
a and b will be omitted from now on unless it causes 
confusion. 

The critical FSS at e = is 



(to) c ~ A e N~ 1/4 , 



(24) 



where A e = -0(0) = (12/6) 1 / 4 r(i)/r(i). We remark 
that e may not be exactly zero at the bulk critical point 
= limjv^oo K c , but may have a finite-size correction 
vanishing exponentially with N. This additional correc- 
tion does not change the leading power-law term in the 
FSS. For e < 0, 



(to) ~ (-eATT 1/2 , 



(25) 



since ip(x) ~ y / —2x/(na) for x — > —00. The scaling 
form in Eq. (|22|) is reproduced from Eq. (|23|) . using the 
limiting behavior of ip(x) ~ ^Jiax/b for x — > 00. The 
crossover between the three scaling regimes occurs at 



and 



-[2/(TcaAl)]N-V 2 , 



[bAl/$a)]N-V\ 



(26) 



(27) 



The scaling behavior of the order parameter (to) is rep- 
resented schematically in Fig. [TJ 

The FSS form for the zero-field susceptibility is given 
by X = N^/^^eN 1 ^; a, b) with 7=1. The scaling func- 
tion <fi(x; a, b) is defined in Eq. (|A5p . 



B. Scale-free networks 

Consider the SF degree distribution P(k) ~ k~ x for 
ko < k < k c and otherwise with the upper cutoff 
k c ~ N 1 ^ and the lower cutoff k = 0(1). We are 
interested in the cutoff exp onent u) > w na t = A — 1 as 
considered in 

H S, where the cutoff-dependent 



<m> ' 


~N" 1/4 














Across Across 



FIG. 1: Schematic plot of {in) versus e in the exponential 
networks. 



FSS in the CP on annealed SF networks. In general, the 
expansion of Eq. (| 14[) for small uq and h is singular as 
z n diverges in the N — > 00 limit for n > A — 1. So one 
should treat the nonanalyticity carefully. Furthermore, 
there is a power-law finite-size correction in the critical 
inverse temperature K c , which plays an intricate role in 
the critical FSS. For all A > 3, the average magnetiza- 
tion is again to ~ a(m) and the magnetic susceptibility 
is identical to that in the exponential networks. 



1 . Finite-size behavior of z n 

As a degree k is an integer, the standard precise ex- 
pression for the degree distribution is 



p(k) = ck- x J2 5 kj, 

j=k 



(28) 



where the normalization factor c is given by c 1 

E-= feo r A with 

z n are given by 



Y2j=ko 3 X with k c = dN 1 /^ . Then, the degree moments 



fee 



j=k 



(29) 



For large k c (large N), we have finite-size corrections 
for the normalization factor as 



- 1 ^c- 1 -fc-( A - 1 )/(A-l), 



(30) 



with c^ 1 = £(A, ko) where the Hurwitz zeta function is 

defined as C(«, I) = E£o0' + P" 

Similarly, we have, up to the leading order in k c , 



,-[A-(n+l)] 

z °° - c ^ 

A - {n + 1) 

Coo fo k c 

,(n+l)-X 



for n < A — 1 

for n = A - 1 (31) 

for n > A — 1 , 



~>+l)-A 
with z%> = CooC(A - n, k ) = ((X - n, fc )/C(A, k ). 
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The critical parameter K c = Z\jz% also has a finite-size 
correction as 

K c ~ [1 + eJ\T a ] , (32) 

with 

K? = z?/z? = C(A - 1, fc )/C(A - 2, fc ), 

e = d-( A - 3 )/[(A-3)C(A-2,fc )], (33) 

and 

a=(A-3)/w. (34) 



2. A > 5 

For A > 5, z n is finite up to n = 4. Hence, the ex- 
pansion of to and /(to) are the same as those in the 
exponential networks up to the order of Uq and up to 
the order of to 4 , respectively, as in Eqs. (fl"8]) and (fT9|) . 
Therefore, their critical behaviors are identical to those 
in the exponential networks, in terms of the parameters 
a, b, e, and N. 

However, unlike the exponential networks, e = (K — 
K c )/K c has a power-law finite-size correction due to the 
TV-dependence of K c . From Eq. (j32|) . one finds that 

e~e b -e f (N) = e b -eN- a , (35) 

where e b = (if — K^°)/K^° is a deviation from the bulk 
critical temperature and a — (A — 3)/w. Therefore, the 
FSS form is given in terms of e b as 

(m(e b , AT)) = N-V*$[(e b - £/)iV 1/P ], (36) 

with /3 = 1/2 and = 2, which shows a simple horizon- 
tal shift of the order parameter curve in Fig. [1] to the 
right (see Fig. [3]). 

The order parameter follows the same scaling law of 
Eq. (|24|) at the ^-dependent pseudo critical temperature 
at e = or e b — ef. On the other hand, at the bulk 
critical temperature at e b = 0, the order parameter is 
given by {fh) b>c = N~ ^^(-eN~ a+1 / i? ). 

For a > 1/2 [uj nat < u < 2(A — 3)], the correction e/ is 
not big enough to shift e^ cro8S = e c T ross + e f with e- ross in 
Eq. (f2"6"|) to cross the bulk critical point e b = [Fig.^a)]. 
Therefore, there is no characteristic change in the critical 
scaling by this shift, except the appearance of a higher- 
order correction to scaling like O (TV^" -1 / 2 ))- 

For a < 1/2 [uj > 2(A — 3)], e^ cross becomes posi- 
tive and both crossovers take place in the side of e b > 
[Fig. [U[b)] ■ The bulk critical point is now in the re- 
gion left to e bcTOSS , where the scaling function behaves as 

4>{x) ~ y/2/(ira) {-xY 1 / 2 . At a = 1/2 [u = 2(A - 3)], 
the scaling variable is finite (x = —efN 1 ^ = — e). 



Therefore, we have the critical FSS at e b = as 

r A.N- 1 ^ (a > 1/2) 

(m) b)C = | A.N- 1 /' 1 (a = 1/2) (37) 

[ B.N-^^/ 2 {a < 1/2), 

where A e = ip{—e) and B e = [2/(7rae)] 1 / 2 with e in 
Eq. (33J. 



3. 3 < A < 5 

For 3 < A < 5, z± and zi are finite, but 2:4 diverges as 
z 4 - k\~ x ~ A/-(5-A)/a> as W ell as 6 = z i (z 1 /z 2 ) 4 . In the 
thermodynamic limit or for uok c + h 3> 1 in finite size 
networks, Eq. (jT5J) has the singular expansion as 

m = h+^uv- ^u*- 2 + 0(h 2 ), (38) 

Zl Zi 

with a constant q ~ cj dxx 1 ~ x (x — tanhx) > 0. For 
uok c + h <C 1 in finite networks, the series expansion 
becomes regular as 

fn = h+—u -—ul + O(ul,h 2 ). (39) 
zi 3zi 

Then, the free-energy density in finite networks is given 

by 

/(m) = -ln2 - ohm - ^em 2 + Q{rh), (40) 



(a) 


<m> 


\ 




-N" 1/4 




. .-1/2 
~(-8N) 












Ecross r "Ef 


Ef Ecross^Ef £|, 


(b) <m>' 


\ 


















N" 1/4 


-1/2 
~(-£N) 









£cross r "Ef Ef Efross^Ef £j, 



FIG. 2: Schematic plot of (m) versus eb in the annealed SF 
networks with A > 5 and (a) oj nat < ui < 2(A — 3) and (b) 
uj > 2(A — 3). Note that the bulk critical point ej, = is 
outside of the critical region in (b). 
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(a) 


<fn> 












•N 4<1+ , 


_ — -"i 


~(-£Np / 
[V 




III 


n 


I 



8* £ 



(b) <m> 




11 ^^S-\ 1 


— 

~£ ra 


~(-£Np 


IV 


; | ~n 2ra £ 2 
III II 


I 


-£5)-£f 


£f £e+£f £*+£f ^ 



FIG. 3: Schematic plots of (m) versus (a) e and (m) versus (b) 
t\, in the annealed SF networks with 3 < A < 5 and uj > w na t- 



where 



Q(m) 



12 



_ 4 

m 



TO*" 1 



for m <C (22/21) k c 1 + h 
for to > (22/21) k^ 1 + h, 



(41) 

with a constant p ~ c/ cfex x (x 2 /2 — In cosh a;) > 0. 
Note that b = {z 1 /z 2 ) 4 z i ~ 6 ^ (5 ~ A)/w with 



&0 = c oo (zrA2°°) 4 rf 5 " A /(5-A) 



(42) 



In the thermodynamic limit, k c becomes infinite and 
the free energy density expansion is singular for all to > 
0. At h = 0, the order parameter scales for e > as 



with 



and 



(to) ~ Ce 13 , (43) 
/9 = l/(A-3), (44) 
C = {z?/z?){z?M\ - \)\Y'^\ (45) 



It is straightforward to calculate the FSS at e = 
by performing the integral in Eq. (fTu| using Eqs. (p!0|) 
and (|4T|) . For cj > w n at, the integral in the region of 
to <C (22/21) fc c 1 dominates and we find 



(m) c ~ 4 s jv"-[i+(5-a)/"]/4 - (67V)" 1 / 4 , (46) 



v / 27(^) (- e )-V2Ar-V2 
, oss ~ -^/(TraA 2 )]^- 1 /^ 



At u> = aj na t , the integrals in both regions contribute, but 
the critical FSS does not change except its amplitude. 

For e < 0, we have (to) ~ 
So, the crossover occurs at e c 
with v- = 2/[l- (5- A)/w]. 

At small positive values of e, we have a nonzero solution 
for to in the region of m <C (22/21 )fc c T 1 as 

(TO) ~ ^JVo AT-(5-^)/(^) e l/2 „ (e/6) V2 ) (4g) 

where e < frN-^-W" with /1 = Coo d-( A - 3 )/[3z^°(5 - 
A)]. For larger e, we have the bulk solution, Eq. (|43|). 
for to in the region of to 3> {zs/zijk^ , where e 3> 
f 2 N-( x - 3 )/ u with / 2 = p(A - l)d-( A - 3 V^f°- 

Hence, there are two crossovers in the side of e > 0. 
The first crossover occurs at 



[b A 2 J(3a))N 



with v +t i = 2/[l — (5 — A)/w], which is the same as is- 
Then second crossover occurs at 



,+.2 



[b C 2 /{ia)]-^/^ N- 1 '^-- 



with i/ +i 2 = w/(A — 3). Note that l/u+ t 2 coincides inci- 
dentally with a in Eq. (135j) . For convenience, we denote 
e ioss b y £ c, eio 2 ss by e», = i>_ by P c , and ^ +j2 by P*. 
They are summarized as 



v c = 2/[l - (5 - \)/u\ 
v* = <jj/(\ — 3). 



(49) 
(50) 



The order parameter is plotted against e in Fig. [3ja), 
where we have one more distinct scaling regime compared 
to the case for A > 5. In (bulk) regime I (e > e„), the bulk 
scaling is valid where the system is free from any finite 
size effect. In (intermediate) regime II (e c < e < e*), 
the system behaves as in a SF network with infinite N 
but with finite k c . In (critical) regime III (|e| < e c ), 
the system feels both finite N and finite k c . Finally, 
the ordinary scaling in the disordered phase appears in 
(disordered) regime IV, where only finite N matters. 

Summing up the results, we need two different scaling 
functions, ip c and ip*, describing the critical region and 
the crossover region to the bulk regime, respectively, for 
the forced cutoff (oj > uJ aa t)- First, near e w 0, we have 



(m(e,N)) = N-P/^ MeN 1 ^"), 



(51) 



A)]. The scaling function 
i (-x)- 1 / 2 for 



with 



with (3 = (lu + 5 - A)/[2(w - 5 
&(0) ~ A s [Eq. (03], Mx) 
x -> -00, and tjj c (x) ~ y/3a/b^ x 1 / 2 [Eq. [jlS]) ] for x -> 
00. Due to the crossover to the bulk regime, this scaling 
function is valid only up to x ~ jVt 1- ^ -1 ^^' 2 , which 
diverges with iV. 

Second, near e w e*, we have 



A a = (i2/6 ) 1 / 4 r(i/2)/r(i/4). 



(47) 



(m(e,N)) =N-' 3 ' P ' ^(e7V 1/p *), 



(52) 
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with (3 = 1/(A — 3). The scaling function ip*(x) ~ 
Cx f3 [Eq. (gSJ)] for a; -> oo and ~ y/3a/b x 1 / 2 

for small a;, but larger than ~ N~^ 1 ~^~ 1 '' u " 2 ; which 
vanishes as AT — ► oo. At w = w nat , the intermedi- 
ate regime II vanishes as e c ~ e* so that the two scal- 
ing functions merge into a single scaling function with 
= j3 = 1/(A - 3) and v c = v* = {\- 1)/(A - 3). 

Now, in terms of the bulk parameter eft, we need to 
replace e by et — e/ in all scaling equations. It implies 
the simple horizontal shift of the order parameter curve 
of Fig. [3] to the right by the amount of e/ = eN~ a 
with a = (A — 3)/<j [Fig. G^b)]. For any to > w aa t, 
we find that e&, CT oss = £ ctoss + e / becomes always posi- 
tive. Therefore, the critical FSS at e& = is (m)b tC = 
N-P/o- ^(-e/iV 1 /^), which results in 



(m)t 



-[l-(A-3)/w]/2 



(53) 



with B e — [2/(7rae)] 1 / 2 . At u = w nat , the critical scaling 
does not change except for its amplitude as (rh)b, c ~ 
2V-V(A-i). 



C. Comparison to quenched networks 

In quenched networks, it is difficult to derive ana- 
lytically the FSS for any model due to the presence of 
quenched disorder. Even in the case that quenched dis- 
order fluctuations are negligible, quenched links gener- 
ate the finite correlations in neighboring nodes, which 
are responsible for the critical point shift (mass shift) 
by a finite amount. This mass renormalization process 
should involve the finite-size correction which determines 
the FSS of the (pseudo-) critical point and the FSS of the 
order parameter follows. 

Recently, Hong et al. [5J conjectured the FSS expo- 
nent based on the droplet-excitation (hyperscaling) ar- 
gument and phenomenological theory. They also numeri- 
cally confirmed that the FSS exponent for the Ising model 
is v — 2 for quenched exponential networks as well as for 
the quenched SF networks with A > 5. For 3 < A < 5, 
v = (A — 1)/(A — 3), regardless of the cutoff exponent u> 
if it is not too strong (u> < A). 

For exponential networks, we find the same FSS for 
annealed and quenched networks. The annealed SF net- 
works with A > 5 exhibit essentially the same FSS as the 
quenched SF networks, but the additional finite-size cor- 
rection on the critical point generates a different FSS 
on the order parameter at the bulk critical point for 
lu > 2(A — 3). For w nat < lu < 2(A — 3), this additional 
correction is irrelevant. 

The annealed SF networks with 3 < A < 5 exhibit 
the anomalous FSS characterized by the combination 
of two different single-parameter scaling functions (or 
two-parameter scaling) with the anomalous intermediate 
regime for any lu > w na t, which is generically distinct 
from the quenched SF networks. However, at to = w nat , 
the intermediate regime disappears and the FSS can be 
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FIG. 4: (Color online) Monte Carlo simulation data at A = 
6. (a) The order parameters (m) c at e = and (m)& >c at 
£b = with different values of (d, u>) = (2.25,5), (3,6), and 
(3, 7) are plotted with symbols. They are compared with the 
analytic results of Eqs. (|24[) and (|37[) which are drawn with 
lines. Numerical values of the coefficients are A e ~ 0.815 860, 
A e ~ 0.687 983, and B e ~ 1.039 20. (b) Scaling piot of 
(m)Af' 3/s versus \e\N 1/i} at ui = 6 with (3 = 1/2 and v = 2. 
Data with different values of TV fall onto the scaling function 
^(|e|A rl/ ' 2 ) drawn with the solid curve. 



described by the ordinary single-parameter scaling func- 
tion with the same exponent v = (A — 1)/(A — 3) as in 
the quenched SF networks. 



IV. NUMERICAL RESULTS 

We performed extensive Monte Carlo (MC) simula- 
tions in the annealed SF networks at various values of 
A and lu to confirm the analytic results in Sec. IIII Bl Es- 
pecially we focus our attention on the cutoff dependent 
FSS behavior. In practice, we consider the Ising model 
on the fully connected network with the heterogeneous 
coupling constants given by the Hamiltonian in Eq. ([8]) 
at all hi = and set J = fcn = 1. 
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FIG. 5: (Color online) Monte Carlo simulation data at A = 4. The parameter values of d are 2.57 for u) — 3 and 1 for the 
other values of u. (a) Numerical data for {m) c and {m}b,c represented with symbols are compared with the analytic results 
represented with straight lines, (b) Scaling plot of (tti)A^ 1 ^ a_1 - ) versus |e|A r( - A_3 - )/ '^ A_1 ' at cu = o; na t = 3. The straight line has 
the slope of f3 — 1/(A — 3) = 1. We also show the scaling collapse of the numerical data at to — 4 according to the FSS forms 
in (c) Eq. (|51[1 and (d) Eq. ([52}. The solid line and the dashed line have a slope of 1/2 and = 1/(A — 3) = 1, respectively. 
We use same symbols of Fig. HJb) and the symbol o for N = 32 768 x 10 3 . 



Using the standard Metropolis single spin update rule, 
we run MC simulations up to 2 x 10 3 - 10 4 MC steps 
for system sizes up to N = 32 768 x 10 3 . The MC data 
are averaged over 100 independent samples of initial spin 
configurations as well as the thermal (temporal) average 
after discarding the data up to 10 3 MC steps for the 
equilibration. 

First, we need to choose a degree sequence {ki} — 
{fci,...,fejv} in accordance with a given degree distri- 
bution P(k) = ck~ x for fc < k < k c as in Eq. l[2"8"j). with 
feo = 0(1) and k c — intfdJV 1 /"], where int[x] denotes 
the integer part of x. Let Nk be the number of nodes 
with degree k. Such a degree sequence can be generated 
deterministically by applying the rule 



N k ' = int 



k'=k 



N £ P(fc') 



L k' = k 



(54) 



to Nk for all k in the descending order from k = k c 



The maximum degree k' c thus obtained may be dif- 
ferent form the target value k c . In fact, k' c can be es- 
timated from the condition NJ2k=k' = 1' which 
yields k' c = dN 1 ^^ + 0(N- 1+U ^ W )] for u > w nat . 
Therefore, Eq. (|5~4"|) indeed yields the degree cutoff scal- 
ing with the prescribed values of d and u) only with a 
higher-order correction. However, when lu = (x> nat , we 
find that k' c = d'N 1 /" with 

d' = d(l + (A - l)d A - x C(A, fco)) _1/(A_1) - (55) 

When one compares numerical data with the analytic 
results, the modified value d' should be used for uj = 
Wnat — A — 1. In this section, we use the degree sequences 
generated deterministically from Eq. ([54")) for various N, 
A, d, and to with fixed fco = 3. 

Monte Carlo simulation data for A = 6 are presented 
in Fig. |U We first test whether the magnetizations (m) c 
at the pseudo critical temperature with e = and (rh)b c 
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at the bulk critical temperature with eb — scale as in 
Eqs. (f2"4"| and (|37[) . respectively. In order to cover the 
three cases of Eq. (|37| . we choose u> = 5, 6, and 7, which 
correspond to a = 3/5, 1/2, and 3/7, respectively. These 
numerical data in Fig. SJa) are in good agreement with 
the analytic results. 

Our analytic theory predicts the full shape of the scal- 
ing function as well as the scaling exponents. We exam- 
ine validity of the FSS form in Eq. (TX2j) in Fig.gjb). We 
present the scaling plot of (rh(e,N)) against ^N 1 /" us- 
ing the Monte Carlo data with A = 6 and uj = 6. These 
data match perfectly well with the analytic curve for the 
scaling function -0(|e|7V 1/2 ) in Eq. (|A"3j) . 

We proceed to the case with A = 4, where the FSS 
behavior is more complicated. We first examine the FSS 
of (rh) c at e = and (rh)b,c at eb — 0. They are predicted 
to follow the power law given in Eqs. (|4"6"|) and (|5"3"|) , re- 
spectively, when u) > uj na _t- When oj — oj na t, the scaling 
is given by the same power law but with modified ampli- 
tudes. Figure [5]^a) presents the plots of (fh) c and (rh}b,c 
against TV at a; = 3, 3.5, 4, and 4.5, which agree well with 
the theoretical curves. 

When uj = cj nat , the FSS is governed with the single 
scaling variable eA^^ 3 ^^" 1 ). In Fig. 0(b), we present 
the scaling plot of (to) V 1 /^ 1 ) against | e | A^- 3 )/^- 1 ) 
at A = 4 and uj = 3. A good data collapse supports that 
FSS form with the single scaling variable. We note that 
the bulk scaling behavior (to) ~ e 13 with j3 = 1/(A — 3) 
sets in only for N > 10 6 . 

We also examine the FSS behavior at uj = 4 (> w nat ). 
Here, the FSS is governed with two scaling variables 
eN 1 ^" and eN 1 ^* . Hence, one cannot expect a data 
collapse over the whole regions in a scaling plot. We first 
test the scaling form of Eq. (f5"Tj) , which is valid in regimes 

II, III, and IV. The scaling plot of {m)N^I Se against the 
scaling variable \e\N x l Va is presented in Fig. [5^c). We 
observe a reasonably good data collapse in regimes II, 

III, and IV except for small network sizes. In the e > 
side, the numerical data align along a straight line of 
slope 1/2, which reflects the scaling (rh) ~ e 1 / 2 in regime 
II. However, they begin to deviate from the straight line 
systematically for N > 4096 x 10 3 as the scaling variable 
increases. This is due to the crossover to regime I. 

Finally, we test the scaling form of Eq. ([52]). which is 
valid in regimes I and II. Figure 0(d) shows the scaling 
plot of (to)V/ 1/ * against the scaling variable \e\N x / v '. 
As expected, we do not have a data collapse for e < 0. 
The data in regimes I and II do not collapse well either. 
The order parameter scales as (rh) ~ e 1 / 2 in regime I 
and then (to) ~ = g 1 i n regime II. Comparing 

the numerical data with the straight lines of slopes of 
1/2 and 1, one finds the signature of the crossover for 
V > 4096x 10 3 . This suggests that the poor data collapse 
may be due to a finite size effect. The system does not 
reach the scaling regime I even at V = 32 768 x 10 3 yet. 



N=1x10^ 




12 3 4 

x = <m> c / [<m> c ] 

FIG. 6: (Color online) The histogram for (m) c /[(m) c ], the or- 
der parameter at e = normalized with the ensemble average, 
where network parameters are A = 4.0, ko = 3, d = 1, and 
uj = ujnut = 3. The curves from different values of V collapse 
onto a single curve. 



V. SAMPLE-TO-SAMPLE FLUCTUATIONS 

In the previous section, we tested the FSS theory for 
the power-law degree distributions generated determin- 
istically from Eq. (|54p . The other way is to draw proba- 
bilistically N values of the degree independently in accor- 
dance with the target distribution function P(k). This 
is adopted in the configuration model [111, 0|. In the 
probabilistic method, the degree sequence varies from 
sample to sample, hence an ensemble average is neces- 
sary. One interesting issue is whether physical quantities 
have the self-averaging property ^20j] against the sample- 
to-sample fluctuations. For finite systems, a sample with 
{hi} = {fcx, . . . , fcjy} drawn probabilistically may show 
the degree distribution P(k) = J2i^k,ki/N, which devi- 
ates from the target distribution function P(k). Then, it 
follows that the degree moments z n = JV kf/N show the 
sample-to-sample fluctuations, purely from the sampling 
disorder. 

Using the same techniques used in our previous publi- 
cation for the CP model [l(| (see Sec. V therein), it is 
straightforward to show that the relative fluctuation R n 
is given by 

d _ [4j - W? 1 /<2n) \ 

Rn = ^^ = N{-W- 1 )> (56) 

where [• ■ ■ ] denotes the sample (disorder) average and 
(n)o = ^2kk n P{k). For exponential networks, all (n) 
are finite, so all degree moments z n are strongly self- 
averaging (R n - V- 1 ) m. 
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In the SF networks with P(k) given in Eq. (l28j) . 

N { n -X+i)/u for n> x-l 
[n}„ \ log A 7 for n = X-l (57) 

0(1) for n< A- 1, 



which leads to 



^-i + (A-i)/ w 



for n > A — 1 



jy-l+(2n-A+l)/a» for (A - l)/2 < n < A - 1 



N~ 



for n < (A - l)/2, 



(58) 

where there are log corrections at n = A — 1 and (A — 1)/2. 
By definition, R n is strongly self-averaging for n < (A — 
l)/2 and is weakly self-averaging for n > (A — 1)/2 except 
that R n is not self-averaging only when uj = w na t for 
ii > A — 1. For example, is not self-averaging for A < 5 
with the natural upper cutoff. 

The relevant quantities involving the degree moments 
are K c — z\/z2, a — zf/z2, and b — (z\j z-i^za,. It 
implies that the critical point location and a are strongly 
self-averaging for A > 5 and at least weakly self-averaging 
for A > 3. However, b is not self-averaging for 3 < A < 5 
with co = Wnat, which determines the amplitude of the 
order parameter in various scaling regimes [see Eqs. (|46[) 
and ([48]) ]. Therefore, we expect widely scattered data 
for the order parameter, depending strongly on sampled 
degree sequences, for 3 < A < 5 with uj — w nat . 

Numerical data are presented to verify the non-self- 
averaging property of the order parameter at e = 0, 
(m) c ~ (bN)- 1 / 4 [see Eq. flU], for the annealed SF net- 
works for A — 4 with to = o; nat = 3. It should not be self- 
averaging because it involves the parameter b. We have 
measured the order parameter (m) c in many samples and 
constructed a histogram of the quantity (m) c /[(m) c ], the 
order parameter normalized with its mean values. 

Figure [6] presents, thus, the obtained histogram. The 
histogram does not sharpen at all, but collapses onto a 
single curve as N increases. This proves the non-self- 
averaging property. 



The FSS exponent is given by v — 2, with which the scal- 
ing variable for the FSS is given by eA 1 /". The scaling 
behaviors in the critical regime [eN 1 '" = 0(1)], in the 
supercritical regime (e N^ 1 ^), and in the subcritical 
regime (e <gC — TV -1 /") are summarized in Figs. Q] and [2] 
For power-law degree distributions with 3 < A < 5, the 
degree cutoff k c ~ N 1 ^ matters and there exist two dis- 
tinct scaling variables eN 1 ^" with D c = 2/[l — (5 — A)/cj] 
and eN 1 /"* with = u>/(\ — 3) when uj > uj na t = A — 1. 
At uj = w na t, the two scaling variables merge into a sin- 
gle one. The scaling behaviors in the supercritical regime 
I (e > A^ 1 / 5 *), the intermediate regime II (Af~ 1/P <= < 
e < N~ x / V *), the critical regime III [eA 1 /^ = 0(1)], and 
the subcritical regime IV (e <C ~N 1 / Vc ) are summarized 
in Fig. [3] The crossover from regime I to II is originated 
from the fmitcncss of the degree cutoff k c , while the crit- 
ical FSS in regime III is from the fmiteness of both k c 
and N. 

The CP on the annealed SF network studied in Refs. [§, 
IToj is also characterized with the two w-dependent FSS 
exponents when 2 < A < 3 and uj > a> n at- The similarity 
between the equilibrium Ising model and the nonequi- 
librium CP suggests that the two-parameter scaling is a 
generic feature of critical phenomena in annealed scale- 
free networks. 

Extensive studies during the last decade have revealed 
that critical phenomena on quenched networks and an- 
nealed networks are characterized with the same set of 
bulk critical exponents such as the order parameter expo- 
nent and susceptibility exponent. However, they display 
distinct FSS behaviors. Annealed networks are charac- 
terized with two FSS exponents, which depend on A and 
lu. In comparison to annealed networks, quenched net- 
works have a quenched disorder in structure. Besides, 
dynamic degrees of freedom on quenched networks have 
finite correlations. It is another big challenge to under- 
stand how these two ingredients cause the distinct FSS 
behaviors, some of which is under investigation plj . 
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We have investigated the FSS of the Ising model on 
annealed networks. The model is mapped to the Ising 
model on a globally connected network with heteroge- 
neous couplings, which allows us to derive the free-energy 
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APPENDIX A: SCALING FUNCTIONS 

From the free energy density function in Eq. (|19p for 
the exponential networks and also the SF networks with 
A > 5, one can easily derive (fh) for small h as 

m ^ ( 12 \ 1/4 U(l/2,r)+a'N^hU(3/A,r) 
\bN J U{l/A,r)+a>N?>/ 4 hU{l/2,ry 
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where a! = (Yljbf^a, r = (3a 2 /bf^eN 1 / 2 , and 

U(s,r)= dy y 3 ' 1 exp(-y + r^/y). 
Jo 

With this, we find the order parameter scaling at h = 

as 

(m(e, AO) = N-^^ieN 1 / 9 ; a, b), (A2) 
where (3 — 1/2, £/ = 2, and 

7/ M ^12V /4 U(l/2,r x) 
^ ;a ' 5)= UJ £7(1/4,^) ' (A3) 

with r = {Zc^/b) 1 ' 2 and x = eN 1 ^. 

The zero-field susceptibility x = {d {rh)) / (dh)\j l={) is 



x{e,N)=W/*4>(fiN 1 !'>-,a y b), 
where 7 = 1 and 



(A4) 



(x; a,b) = a 
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1/2 



U(3/A,r x) U 2 (l/2,r x) 
U(l/4,r x) [/ 2 (l/4,r x) 



Using the properties of the function U(s,r) such as 



U(s,r) 



2F(2s) 

r( s ) 



-2.s 



1 + 



4(s-l)(s-2) 



(r — > — oo) 
(r = 0) 
(r -> oo), 



(A5) 



(A6) 

one can show tp(x) ~ ^/3a/6 x 1 / 2 or ^/2/ (7ra)(— .x) 1/>2 
for x — > ±oo, and </>(x) ~ (2x)^ 1 or (1 — 2/ir)(— x) _1 
for x — > ±oo. We remark that the usual magnetization 
and the magnetic susceptibility become m ~ a(m) and 
X - ax- 



For the SF networks with 3 < A < 5, the scaling func- 
tion i/j c (x) for the order parameter near e w in Eq. (|5T|) 
behaves in the same way as the above ip( x ) except for 
replacing j3 by 0, v by 9 C , and b by 6o- For example, r 
becomes r = (3a 2 /6 ) 1/2 eA^ 1/Pc . The susceptibility scal- 
ing function also changes in the same way. 
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